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Abstract 

We explore the phenomenology of microscopic black holes in the S 1 fZz Randall-Sundrum (RS) 
model. We consider the canonical framework in which both gauge and matter fields are confined 
to the brane and only gravity spills into the extra dimension. The model is characterized by two 
parameters, the mass of the first massive graviton (mi), and the curvature l/£ of the RS anti-de 
Sitter space. We compute the sensitivities of present and future cosmic ray experiments to black 
hole mediated events, for a wide range of £ and m\, and compare them with the sensitivities of 
Tevatron Runs I and II to higher-dimensional physics. As part of our phenomenological analysis, 
we examine constraints placed on £ by AdS/CFT considerations. 
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One of the most exciting predictions of sub-millimeter extra dimensions |T], |2| is the pro- 
duction of black holes (BHs) in particle collisions with center-of-mass energy larger than a 
TeV and sufficiently small impact parameter [|3], £|, |5], |6|, 0, ||. Although colliders have not 
yet attained the energies required to probe this new strong quantum gravitational effect, 
the extraordinarily high center-of-mass energies achieved at the top of the atmosphere in 
ultra-high energy cosmic ray collisions are high enough to render any change in spacetime 
dimensionality detectable [TU], |TT], [I2|, [T3|. Among these cosmic rays, a nearly guaranteed 



flux of neutrinos (produced through interactions of extremely high energy protons with the 
cosmic microwave background) could produce BHs which decay promptly initiating deeply 
developing air showers far above the Standard Model (SM) rate |J, and with very distinctive 
characteristics [ 1Q| . In addition, neutrinos that traverse the atmosphere may produce BHs 



through interactions in the ice or water and be detected by neutrino telescopes [13| . More 



over, in scenarios with asymmetric compactifications the production of brane configurations 
wrapped around small extra dimensions may be competitive with BH production |14| . Very 
recently, based on the absence of a significant signal of deeply developing showers reported 
by the AGASA Collaboration we derived new limits on the fundamental Planck scale |12] . 



In this paper we expand upon this study and examine in more detail the phenomenological 
implications of BH production in the Randall-Sundrum (RS) scenario 0. 

The RS model consists of two 3-branes (with equal and opposite tensions cr pianck = —u SM = 
12M 3 /£) which rigidly reside at S 1 /Z% orbifold fixed points at the boundaries (y = and 
y = 7rr c ) of a slab of anti-de Sitter (AdS) space of radius I. The classical action describing 
the system is given by JJ] 

S = M 3 [ d A x r° dy^=g + + I d A x a SM + / d 4 x yf^q a planck , (1) 



where R is the 5-dimensional Ricci scalar in terms of the metric g^ U) M is the fundamental 
scale of gravity, and rjij is the flat Minkowskian metric. In what follows, the Latin subscripts 
extend over ordinary 4-dimensional spacetime, whereas Greek subscripts represent all 5 
dimensions. The metric satisfying this Ansatz (in horospherical coordinates) reads 

ds 2 = e~ 2lvl/e rj i:j dx i dx j + dy 2 . (2) 

Examination of the action in the 4-dimensional effective theory leads to || 

Mpj = M 3 t (l-e- 27Trc/e ) , (3) 

where Mp\ is the reduced effective 4-dimensional Planck scale. Now, assuming that SM fields 
are localized on the 3-brane at y = nr c , one finds that a field with the fundamental mass 
parameter m will appear to have the physical mass m = e~ 7TTc ^ e m . Hence, TeV scales can 
be generated from fundamental scales of order Mp\ through the exponential warping factor. 
Specifically, the observed hierarchy between the gravitational and electroweak mass scales 
is reproduced provided r c /£ w 12. The 4-dimensional phenomenology of this model (only 
gravity propagates into the bulk) is governed by two parameters: c = (£Mpi) _1 , which is 
expected to be near though somewhat less than unity, and m% which is the mass of the first 
Kaluza-Klein graviton excitation [15[ . 



Two different types of BHs can be produced in trans-Planckian particle collisions within 
this set up: (i) AdS/Schwarzschild BHs that propagate freely into the bulk (generally falling 
towards the AdS horizon once produced) and (ii) tubular pancake shape BHs that are bound 
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to the brane 0, |16 . To study the phenomenology of the latter, it is convenient to define a 
new variable z = £e y ^. In such a coordinate system the metric 
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ds 2 = — (dz 2 + rjij dx 1 dx^ (4) 

is conformal to a 5-dimensional flat metric. To describe the SM brane we introduce the 
coordinate w = z — z c , where \w\ E (0, w c ), w c = £ [e KTc l l — 1), and the TeV brane is located 
at w = 0. After a conformal redefinition 




one obtains W? 



9^ = k— ^ I ft- («) 



*=i^V*-£. («) 



where i? is the Ricci scalar calculated with the metric g^ v . Initially Minkowskian, the metric 
g will be modified, under conditions to be delineated, to include as a patch the 5-dimensional 
Schwarzschild solution. The gravity sector can now be rewritten as 



^gravity = / d^X \ dw J -Q | | | R 



where 



o \z c + wj \ (z c + w) 2 J 

M 3 I d 4 x l o dw^f-g I ^-^2 J > ( 7 ) 



M = Mz c /£ . (8) 



It is worth emphasizing that M refers to the canonical frame, and is expected to be of order 
1 TeV. 

We will be interested in the domain of parameter space for which a high energy collision 
as viewed from the SM brane can result in the formation of a 5-dimensional spherical flat 
space black hole. A calculation of metric perturbations due to a source on the w = brane 
can be made in the flat 5-dimensional space-time approximation, i.e., ignoring the effects of 
the AdS term in Eq. (^) , if w z c , where z c is the AdS curvature as viewed on the SM 
brane, with a gravitational constant M. This implies that to use flat space BH formulae we 
must require r s <C z c , where r s denotes the size of 5-dimensional Schwarzschild radius in the 
canonical frame [1^ 



bhj 



2Mbh 

^SvrM 3 ' [) 

where Mbh = V§ is the BH mass, s is the center-of-mass energy in terms of g^ u , and 
M D = (4tt) 1 / 3 M. 1 From Eqs. ©, (@), and the forg oing relation between Mjj and M we 



1 We set = (4Gs) 1 , where G5 is the 5-dimensional Newton constant. 
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find that the condition r s < z c leads to an upper bound on the mass of the black hole for 



which this picture is valid [20 



M BH /M D < 24 c 



-4/3 



(10) 



where, once more, c = (iM pi) . When the energy exceeds this bound, the behavior of the 



cross section may be analyzed within the AdS / CFT dual picture |21|] , and may assume the 
ln 2 E behavior conforming to the Froissart bound p2]j . 

As mentioned earlier, the parameter c is expected to be small. Comparison of the RS 
brane tension and the D-brane tension in perturbative heterotic string theory suggests that 
c < 0.1 [Tj|. Work in recent years linking the Randall- Sundrum brane- world mechanism to 
the conjectured AdS/CFT correspondence allows an independent and similar bound on 
c. The one-sided brane world relation, Eq.(BI) with r c — > oo, is equivalent to 



{Ml) 



-3/2 



(11) 



In this "single brane" limit f23fl , Duff and Liu |24| pointed out a complementarity be- 
tween AdSs and an Af = 4 superconformal U(AQ Yang-Mills theory living on the four- 
dimensional brane (boundary): corrections to Newton's law from both sectors are equal 



when the AdS/CFT Weyl anomaly relation [25 

2G 5 iV 2 = 



(12) 



or 



(Ml) 3 = N 2 /(8n 2 ) (13) 

holds. For finite but large r c , deviations from conformality are exponentially damped in 
the infrared ||26|| . BH creation in SM brane collisions becomes significant if a weak-gravity 
5-dimensional description is valid [57]. This occurs for Ml 3> 1, and hence c < 1. A more 
quantitative estimate emerges by first combining Eqs 



TTD and (O) to obtain 



c = 2V2n/N 



(14) 



An approximate lower bound on N may be surmised by noting that the KK gravitons find a 
dual description as glueballs resulting from strong coupling in the (approximate) CFT sector 
||27|| . In the large- N analysis of 't Hooft P8f , this requires that the planar loop expansion 
parameter gY M N/8n 2 > 1. For (yfyivi/^ 71 " — 0.1, we obtain from ( |i~4"P 



< 0.1 



(15) 



The key question now is at what mass ratio (Mbh/^Td) is the BH description valid. 
According to the semiclassical prescription, the BH evaporation is governed by its Hawking 
temperature 

Th = ^- (16) 



2vrf s 

Since the wavelength A = 2n /Th corresponding to this temperature is larger than the BH 
size, to a very good approximation the BH behaves like a point-radiator with entropy 



S = - n M BB r s 



/32tt /M] 
~27 \~M 



3/2 



1311 



(17) 



D 
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The magnitude of the entropy indicates the validity of this picture. Thermal fluctuations 



due to particle emission are small when S 1 [29], and statistical fluctuations in the micro- 
canonical ensemble are small for 3> 1 ||. In searches for BH mediated events at colliders, 
it is essential to set x m i n = M^/Mjj high enough that the decay branching ratios predicted 
by the semiclassical picture of BH evaporation are reliable. The QCD background is large, 
and therefore the extraction of signal from background at hadron colliders depends on know- 
ing the BH decay branching ratios reliably. This is especially true if one is attempting to 
determine discovery limits, where the overall rates for BH production are not necessarily 
large. Thus, in collider searches, a cutoff of s min = 5.5 (i.e., S > 25) or more may be ap- 
propriate. By contrast, the search for deeply penetrating quasi-horizontal showers initiated 
by BH decays can afford to be much less concerned with the details of the final state, since 
the background is almost nonexistent. As a result, the signal relies only on the existence 
of visible decay products, which, in this context, includes all particles other than neutrinos, 
muons, and gravitons. Indeed, there is very little about the final state, other than its total 
energy and to some degree its multiplicity and electromagnetic component |I0| , that we can 
reasonably expect to observe, since detailed reconstruction of the primary BH decay process 
is not possible at cosmic ray detectors. It seems reasonable to choose a significantly lower 
value of Mgg 1 than is needed for collider searches; in our estimates of rates for cosmic ray 
facilities we will take x m i n as low as 1, or S as low as 2. While BHs of mass around Md will 
be outside the semiclassical regime, it seems quite reasonable to expect that they will nev- 
ertheless decay visibly, whatever 5-dimensional quantum gravitational description applies. 
Finally, as a consequence of Eq.(|T0"D there are upper bounds on the entropy. Combining the 
latter with Eq. (|17|) we obtain 

S < 250 cT 2 . (18) 

To model the collision between two partons i and j, with center-of-mass energy a/J and 
impact parameter b ^ f s , we consider a scattering amplitude of an absorptive black disk 
with area tc f 2 s || . Criticisms |30| of the assumptions leading to these cross sections have been 



addressed |3T| . Additionally, when considering non-zero impact parameters in 5 dimensions, 



one must be sure that the production of BHs is favored over that of black rings [[32] for the 
relevant values of the angular momentum J. For a 5-dimensional Kerr solution [19| with a 
single angular momentum J (taken in the direction of the brane 0), the horizon size and 
J for a mass Mbh are related: 



4 2 _, J 32G 5 M BH , |() . 
9 rfc+ M| H " 27vr • 



A geometric cross section requires absorption in angular momenta up to J ~ On 
incorporating this condition in (|l9|) one finds that the maximum value of J required for the 



geometric cross section is 0.83 Mbh a/32 G5 Mbh /27ir. This is below the lower limit on the 



spin of a black ring, 0.92 Mbh y 32 G 5 M B h/277t |3"2fl . Therefore, it is consistent to consider 
only the production of 5-dimensional BHs with geometric cross section. For the same reasons 
as discussed in |12| , we will calculate our cross sections with the Schwarzschild radius. 

In our investigation of BH production by cosmic rays, we will be most interested in 
collisions of neutrinos with atmospheric nucleons. In order to obtain the vN cross section 
we take the geometric cross section a, fold in the appropriate parton densities and integrate 
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FIG. 1: Allowed region of the RS parameter space. Resonant production of the first KK graviton 
excitation in the Drell-Yan and dijet channels at the Tevatron [^3| has excluded the region to the 
left of the solid bumpy curve |l5|] , whereas an analysis of the oblique parameters excludes the region 
below the smoothly falling curve (dotted) ^4|. The diagonal dashed lines correspond to future 95% 
CL parameter exclusion region at Run II at the Tevatron under the assumption that no signal is 
found |l5|] . The smoothly (dashed-dotted) rising curves indicate contours of constant Mb in the 
c — mi plane. The region above AGASA's curve is excluded at 95% CL. For an observation time 
of 5 (3) years, Auger (OWL) will probe the region above its curve. The unshaded area indicates 
the region delineated by Eq. (|i~5|) . 

over the parton energies || 

a{uN - BH) = V [I dx Oi(v^i) fi(x, Q) , (20) 

where s = 2m^E v ^ the sum is over all partons in the nucleon^and the /, are parton 
distribution functions. We set the momentum transfer Q = min{MBH, 10 TeV}, where the 
upper limit is from the CTEQ5M1 distribution functions [Q. Also, as a consequence of 
Eq. ( [10|) , there will be an upper bound on the parton subenergies. However, the effect is 
negligible over virtually the entire region of c. 
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With this in mind, for an incoming neutrino flux dQ/dE u , the event rate for deep showers 

/d^> 
dE v N A — - a(uN -> BH) A{E V ) T , (21) 
dhj u 



where JS a = 6.022 x 10 23 is Avogadro's number, A(E U ) is the acceptance for quasi-horizontal 
showers in cm 3 water-equivalent steradians, and T is the experiment running time. In what 



follows we adopt the cosmogenic neutrino flux estimates of Protheroe and Johnson |34 
assuming an energy cutoff in the injection spectrum at E cuto s = 3 x 10 21 eV. This is in very 
good agreement with the most recent evaluation of the cosmogenic neutrino flux assuming 
a strong cosmo logical evolution on the cosmic ray sources, which scales like oc (1 + z) 4 for 
redshift z < 1.9 and becomes flat at higher redshifts [[£J. To explore possible additional con- 
tributions from semi-local nucleon sources, we also consider below the cosmogenic neutrino 
flux estimates of Hill and Schramm [36]. Bounds on Mp emerge as a result of adopting the 
AGASA limit of 1 possible event observed, with a background of 1.72 events |ft7| . This places 
a limit of < 3.5 events at 95% CL ||38|| . Inserting our cross section into Eq. (|2T| ) and requiring 
M < 3.5 leads to exclusion of the region M D < 0.70 ± 0.05 TeV at 95% CL |2f . The varia- 
tion in this bound reflects uncertainty in the neutrino flux. The forthcoming facilities of the 



Auger Observatory |39| will probe fundamental Planck scales, Md ~ 1-55 ± 0.15 TeV |T2 



assuming 5 years of data and no excess above SM neutrino background. The variation in 
this case includes uncertainties in the hadronic background and the experimental aperture. 

Beyond Auger, NASA has authorized studies for the satellite cosmic ray detection facility 
known as the Orbiting Wide-angle Light-collectors (OWL) jHJ, projected for 2007. In a 
comprehensive recent study, BH event rates at OWL for n > 2 large extra dimensions 
were reported [^TJ. The greatly increased aperture of this space-based facility will allow the 
detection ~ 10 BH/yr for x m i n as large as 5. This is well into the region of large entropy, 
permitting comparison with characteristic shower profiles. In the case of present interest 
(n — 1), the "eyes of the OWL" will substantially extend the region of the c — mi plane 
probed in cosmic ray observations. Under the assumption that the cosmic ray sources have 
a strong cosmological evolution, the expected background from SM neutrino interactions 
(all flavors) is 3 events/year jy|. On the basis of zero hadronic background, three years of 
observation implies a deviation from the SM at a 95% CL for 7.77 events observed above 
background [j3~8fl - If these are BHs produced with x min = 1, this represents a sensitivity to 
M D < 2.8 TeV at the 95% CL pi. 



Summary. The Randall- Sundrum model predicts a series of TeV-scale graviton resonances 
with weak scale couplings to SM fields: this makes the Tevatron an outstanding probe of RS 
physics. Analysis ]l5j of Tevatron data |43| for anomalous Drell-Yan and dijet production 



as well as the calculation of indirect contributions to electroweak observables |44] already 



place significant constraints on the c — m\ plane. In this work we have developed the 
formalism for calculating RS black hole production through cosmic ray neutrino collisions 
in the atmosphere, and have delineated additional constraints placed on the RS model by 
existing and future cosmic ray data. The complete allowed region is summarized in Fig. 1. 
The lines of constant M D reflect the relation m 1 = 3.83 c 2/3 (4n)' 1/3 M D . AGASA is able 
only marginally to probe the parameter space allowed by the Tevatron in the limited region 
mi ~ 1 TeV and c > 0.6 — 0.7. Five years of running at Auger can reveal events in the 
region mi > 1 TeV and c > 0.5 which will not be tested by the Tevatron Run II A. A year of 
observation by OWL will probe a substantially larger part of the parameter space for small 
c (see Fig. 1). We have also delineated in our figure the range of the parameter c (related 
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to the AdS radius) favored through AdS/CFT considerations. 
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